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Abstract: We study the phenomenology of the CMSSM/mSUGRA with non-thermal 
neutralino dark matter. Besides the standard parameters of the CMSSM we include the 
reheating temperature as an extra parameter. Imposing radiative electroweak symmetry 
breaking with a Higgs mass around 125 GeV and no dark matter overproduction, we 
contrast the scenario with different experimental bounds from colliders (LEP, LHC), cosmic 
microwave background (Planck), direct (LUX, XENONIOO, CDMS, IceCube) and indirect 
(Fermi) dark matter searches. The allowed parameter space is characterised by a Higgsino- 
like LSP with a mass around 300 GeV. The observed dark matter abundance can be 
saturated for reheating temperatures around 2 GeV while larger temperatures require extra 
non-neutralino dark matter candidates and extend the allowed parameter space. Sfermion 
and gluino masses are in the few TeV region. These scenarios can be achieved in string 
models of sequestered supersymmetry breaking which avoid cosmological moduli problems 
and are compatible with gauge coupling unification. Astrophysics and particle physics 
experiments will fully investigate this non-thermal scenario in the near future. 
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1 Motivation for non-thermal dark matter 

One of the main particle physics candidates for dark matter (DM) is a stable neutralino y 
which emerges as the lightest supersymmetric particle (LSP) in several scenarios beyond 
the Standard Model (SM). The DM relic abundance is generically assumed to be produced 
thermally by the following process: the LSP is in a thermal bath in the early universe, 
subsequently drops out of thermal equilibrium and freezes-out at temperatures of order 
Tf ^ 171^120 when DM annihilation becomes inefficient. 

However, we have no direct observational evidence of the history of the universe before 
Big Bang Nucleosynthesis (BBN) for temperatures above Tbbn — 3 MeV. There is therefore 
no reason to assume a very simple cosmological history characterised by just a single period 
of radiation dominance from the end of inflation until BBN. In fact, the presence of a period 
of matter domination between the end of inflation and BBN could completely change the 
final prediction for the DM relic density if the reheating temperature at the end of this 
period of matter dominance is below Tf [1, 2]. 

This non-thermal picture emerges generically in UV theories like string theory due 
to the ubiquitous presence of gravitationally coupled scalars [3-5]. During inflation these 
fields, called moduli, get a displacement from their minimum that is in general of order 
Mp [6]. After the end of inflation, when the Hubble constant reaches their mass, H ^ mmod? 
they start oscillating around their minimum and store energy. Redshifting as matter, 
they quickly dominate the energy density of the universe which gets reheated when the 
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moduli decay. Being only gravitationally coupled, the moduli tend to decay very late when 
H ^ T ^ '^mod/^p* corresponding reheating temperature 


Tr ^ ^/TMp ^ rumod 



( 1 . 1 ) 


has to be larger than Tbbn in order to preserve the successful BBN predictions.^ This 
requirement sets a lower bound on the moduli masses of order rumod ^ 30 TeV [3]. 

Generically in string compactifications supersymmetry (SUSY) breaking effects de¬ 
velop a mass for the moduli and generate by gravity or anomaly mediation soft terms 
of order Mgoft. Due to their common origin, the mass of the lightest modulus mmod is 
therefore related to the scale of the soft terms as Mgoft = ^n^mod- Given the cosmologi¬ 
cal constraint rumod ^ 30 TeV, only models with p <C 1 can allow for low-energy SUSY 
to solve the hierarchy problem. Values of p ^ 0(10“^) can come from loop suppression 
factors [8-10] while much smaller values k ^ 0(10“^ — 10“^) can arise due to sequestering 
effects [11, 12]. For Mgoft ^ C^(l) TeV, the corresponding reheating temperature becomes 
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which for 10“^ ^ ^ ^ 10“^ is between 0(10) MeV and 0(10) GeV. This is below the 
freeze-out temperature for LSP masses between 0(100) GeV and 0(1) TeV which is Tf ^ 
0(10 — 100) GeV. Therefore any DM relic density previously produced via the standard 


thermal mechanism gets erased by the late-time decay of the lightest modulus. In this new 


scenario, the LSP gets produced non-thermally from the modulus decay. 

From a bottom-up perspective, non-thermal cosmological histories can also enlarge the 
available parameter space of different DM models consistent with direct and indirect de¬ 
tection experiments, due to the presence of the additional parameter Tr. This is appealing 
as it is very hard to reproduce a correct thermal relic density in the CMSSM/mSUGRA 
(see for instance [13]) since a Bino-like LSP tends to overproduce DM (apart from some 
fine-tuned cases like stau co-annihilation and A-funnel or in the case of precision gauge 
coupling unification [14]) while for a Higgsino- or Wino-like LSP the relic density is in 
general underabundant (except for cases like well tempered Bino/Higgsino or Bino/Wino 
DM [15]). 

The purpose of this paper is to study non-thermal DM in the CMSSM/mSUGRA 
where the free parameters are: the standard parameters of the CMSSM/mSUGRA [16], 
i.e. the universal scalar mass m, gaugino mass M and trilinear coupling A defined at the 
GUT scale, tan /3 and the sign of /i, with in addition the reheating temperature Tr from the 
decay of the lightest modulus. We shall follow the RG running of these parameters from 
the GUT to the electroweak (EW) scale and require a Higgs mass rrih — 125 GeV, a correct 
radiative EW symmetry breaking (REWSB) and no DM non-thermal overproduction. We 
shall then focus on the points satisfying all these requirements and we will impose on them 
several phenomenological constraints coming from LEP [17], LHC [18], Planck [19], Fermi 

^Tr has also to be lower than the temperature above which the internal space decompactifies [7]. 
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(pass 8 limit) [20], XENONIOO [21], CDMS [22], IceCube [23] and LUX [24]. Moreover 
we shall focus only on cases where the LSP has a non-negligible Higgsino component since 
Bino-like DM requires a very low reheating temperature which is strongly disfavoured by 
dark radiation bounds in the context of many string models [25]. Interestingly we shall 
find that the constraints from Fermi and LUX are very severe and do not rule out the 
entire non-thermal CMSSM parameter space only for reheating temperatures Tr > 0(1) 
GeV. The best case scenario is realised for Tr — 2 GeV where a Higgsino-like LSP with a 
mass around 300 GeV can saturate the observed DM relic abundance. For larger reheating 
temperatures the LSP Bino component has to increase, resulting in strong direct detection 
bounds which allow only for cases with DM underproduction. Values of Tr above 1 GeV 
require values of ^ 0(10“^ — 10“^) which can be realised only in models where the 
CMSSM is sequestered from the sources of SUSY breaking [11, 12]. Apart from DM, these 
models are very promising since they can be embedded in globally consistent Calabi-Yau 
compactifications [26], allow for TeV-scale SUSY and successful inflationary models [27], 
do not feature any cosmological moduli problem,^ are compatible with gauge coupling 
unification and do not suffer from any moduli-induced gravitino problem [29] . 

In Sec. 2 we discuss CMSSM soft terms, in Sec. 3 we analyse the non-thermal CMSSM, 
in Sec. 4 we discuss our results and conclusions are given in Sec. 5. 

2 CMSSM soft terms 

In a UV completion of the MSSM like string theory, SUSY is spontaneously broken by 
some dynamical mechanism which generates particular relations between the soft terms 
via gravity, anomaly or gauge mediation. In the case when the soft terms are universal at 
the GUT scale, they are given by the scalar mass m, the gaugino mass M, the trilinear 
coupling A and the bilinear Higgs mixing B. We can generically parameterise these soft 
terms and the fi parameter at the GUT scale as: 


m = a |M| , A — bM^ B = cM^ jj. = d M ^ 


( 2 . 1 ) 


where, in a stringy embedding, the coefficients a, 5, c and d are functions of the underlying 


parameters while the gaugino mass M sets the overall scale of the soft terms in terms of 


the gravitino mass m^/2' order to perform a phenomenological analysis of this scenario 
one has to follow the renormalisation group (RG) evolution of these soft terms from the 
GUT to the EW scale and impose the following constraints: a correct REWSB, a Higgs 
mass of order rrih — 125 GeV, no DM overproduction and no contradiction with flavour 
observables and with any experimental result in either particle physics or cosmology. 

A viable REWSB can be obtained if at the EW scale the following two relations are 


satisfied: 



( 2 . 2 ) 


^Ref. [28] provides a significant bound on moduli masses and the number of e-foldings during inflation 
which can be a challenge for many models. 
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where: 


(2.3) 


sin(2^) = 


2\Bi^\ 


mi 


+ m% 

Blu 


Given that the requirement of a correct REWSB fixes only the magnitude of ji leaving its 
sign as free, the parameters (2.1) are typically traded for the standard CMSSM/mSUGRA 
parameters: 

{H^u) 


m — a \M \, A — b M ^ tan /3 = 


(^5) ’ 


sign(M), 


(2.4) 


where one runs m, M and A (or a, b and M in our case) from the GUT to the EW scale 
with a particular choice of tan/3 and sign(/i). Eqs. (2.2) and (2.3) then give the value of 
B and fi at the EW scale. This is the way in which typical spectrum generators operate.^ 
The boundary values of B and /i at the GUT scale which give a correct REWSB can be 
obtained by running back B and fi from the EW to the GUT scale. In this way we obtain 
the values of the coefficients c and d. In a viable UV model, these values of c and d have 
to be compatible with the values allowed by the stringy dynamics responsible for SUSY 
breaking and the generation of soft terms. 


3 Non-thermal CMSSM 


As motivated in Sec. 1, we shall consider scenarios where the LSP is produced non-thermally 
like in the case of string compactifications where the reheating temperature Tr from the de¬ 
cay of the lightest modulus is generically below the thermal freeze-out temperature [4, 5]. 
This reheating temperature represents an additional parameter which has to be supple¬ 
mented to the standard free parameters of the CMSSM (a, 6, M, tan/3 and the sign of 
/i). We call this new scenario the ‘non-thermal CMSSM’ which is characterised by the 
following free parameters: Tr, a, 6, the gaugino mass M, tan/3 and the sign of /i. 


3.1 Non-thermal dark matter relic density 

The abundance of DM particles y produced non-thermally by the decay of the lightest 
modulus is given by [2] : 


/^\NT^ (aann^)P [g^jTf) fTf\ 

V W [U Job. {a,r.nV)f V 9*iTR) \TrJ 




(3.1) 


where is the number of relativistic degrees of freedom, (o’ann'^) ^ 2 x 10“^^cm^ s“^ is 

the annihilation rate at the time of freeze-out needed in the thermal case to reproduce the 
observed DM abundance: 



^Here we use SPheno v3.3.3 [47]. 
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whereas the yield of particle abundance from modulus decay is: 


■5^od — 


STr 

4777/jnod 


^mod 

Mp 


(3.3) 


Br^ denotes the branching ratio for modulus decays into i?-parity odd particles which 
subsequently decay to DM. 

The expression (3.1) leads to two scenarios for non-thermal DM: 


1. ‘Annihilation scenario’: in this case the DM abundance is given by the first term on 
the right-hand side of eq. (3.1) and the DM particles undergo some annihilation after 
their initial production by modulus decay. In order to avoid DM overabundance one 
needs 

(«)/>(w>j'’y^i^(|i). (3.4) 

Given that Tr < Tf and g^{TR) < g^{Tf)^ this scenario requires (o’ann'^)/ > (^^ann'^^) 
as in the case of thermal underproduction. This condition is satisfied by a Higgsino- or 
Wino-like LSP but not by a pure Bino-like LSP which would generically lead to non- 
thermal overproduction (apart from the aforementioned cases). However, given that 
we shall focus on models with just gravity mediated SUSY breaking (contributions 
from anomaly mediation are subleading) and universal gaugino masses at the GUT 
scale as in [11, 12], the LSP can never be Wino-like due to the RG running of the 
gauginos.^ In this context, the ‘Annihilation scenario’ requires a Higgsino-like DM. 
Let us finally point out that the non-thermal DM relic density can be written in 
terms of the thermal one as: 


qNT, 2 _ 1 9*{Tf) f R \ 2 .3 54 

For 5 GeV < Tj < 80 GeV (corresponding to 100 GeV < < 1.6 TeV), the top, 

the Higgs, the Z and the are not relativistic, giving g*{Tf) = 86.25 and: 






0.142 


I 10.75 

9*{Tr) 





(3.6) 


2. ‘Branching scenario’: in this case the DM abundance is given by the second term 
on the right-hand side of eq. (3.1) and the DM particles are produced directly from 
the modulus decay since their residual annihilation is inefficient. In this case both 
large and small cross sections can satisfy the DM content since the annihilation cross 
section does not have any impact on the DM content. This scenario is very effective 
to understand the DM and baryon abundance coincidence problem [31]. Given that 
in general we have Br^ > 10“^, in order not to overproduce DM for LSP masses of 

^Moreover a Wino-like LSP has a significantly larger annihilation cross section than a Higgsino-like LSP 
resulting in a strong conflict with Fermi bounds for sub-TeV Wino-like DM [30]. 
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order hundreds of GeV, one needs l^mod ^ 10 This condition requires a very low 
reheating temperature: 

Tr < 10-V„,od = . (3.7) 

For Mgoft ^ 0(1) TeV and k. ^ O(10~^ — 10“^) we find Tr < (9(10) MeV. In order 
to obtain such a low reheating temperature one has in general to consider models 
where the modulus coupling to visible sector fields is loop suppressed [5]. However 
in this case it is very challenging to avoid a large modulus branching ratio into 
hidden sector light fields like stringy axions [32] and so typically dark radiation is 
overproduced [25]. Therefore the ‘Branching scenario’ does not seem very promising 
from the phenomenological point of view. 

3.2 Collider and CMB constraints 

Due to the considerations mentioned above, if the LSP is Bino-like we generically get DM 
overproduction also in the non-thermal case. We shall therefore look for particular regions 
in the non-thermal CMSSM parameter space where the LSP has a non-negligible Higgsino 
component. We have developed a Monte Carlo programme to find the regions of this 
parameter space where the LSP is Higgsino-like, the Higgs mass is around 125.5 GeV,^ 
REWSB takes place correctly and the following phenomenological constraints are satisfied: 

• LEP [17] and LHC [18] constraints on neutralino and chargino direct production: 
mx ^ 100 GeV; 

• LHC [18] bounds on gluino and squark masses: rrig > 1300 GeV; 

• LHC constraints from flavour physics: M^M~) [33] and the constraint on 

BR(5^ 57 ) [34]; 

• Planck data on DM relic density [19]. 

To avoid complications with the applicability of the standard Spheno version, we restrict 
our scan on the following parameter ranges: tan/3 = 1 to 55, a = OtolO, 5 = —5to5, 
and the universal gaugino mass at the high scale M — 0.3 — 3 TeV. The results are shown 
in Figs. 1-3 for positive /i (the LSP relic density has been calculated using micrOMEGAS 
v3 [48]). 

The plots in Fig. 1 show the points surviving the above constraints in the AjM-vnlM 
plane (at the GUT scale). The points fit into a V-shaped band illustrating a slight hierarchy 
between scalar and gaugino masses (m > 5M) and values of A almost symmetric around 
A ^ M. The regions shown in the plot are mostly for Tr rather smaller than Tf which 
keeps mostly the focus point regions in the allowed parameter space. The coannihilation 
and A-funnel regions can also contribute to the allowed parameter space but they are very 

^Both ATLAS and CMS give values of the Higgs mass between 125 and 126 GeV. In what follows, we 
will consider ranges of values in this region, allowing to some extent for the uncertainty in the spectrum 
generators as well. Allowing for a larger uncertainty in the Higgs mass does not alter the following results 
qualitatively. 
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Figure 1. Correlation between a = m/M and h = A/M ioi different LSP compositions (left) and 
DM thermal relic densities (right) in the region where the LSP is at least 10% Higgsino. 


fine tuned. The V-shape of our plots is caused by the focus point region which can be 
obtained by setting fi ^ mz in the EWSB condition with loop corrections. In fact, the 
dependence of on M and m arises through which depends on the UV soft 

terms in the following way: + g[ff)AjM + h{Q){AlM)‘^ + e(Q)(m/M)^), where 

/, h and e depend on dimensionless gauge and Yukawa couplings (e also includes the 
tadpole correction from the stop loop) and Q is the SUSY breaking scale [35]. A leading 
order cancellation in this expression, as needed to achieve a small /i-term in (2.2), gives a 
V-shaped band in the AjM-mjM plane. We also apply the Higgs mass constraint in this 
parameter space which depends on the square of Xt = At — g cot /3 [36] and preserves 
the V-shape due to its dependence on [A/M^. 
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Figure 2. Correlation between a = m/M and b = A/M for different values of tan/3 (left) and 
gaugino mass (right) in the region where the LSP is at least 10% Higgsino. 
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To illustrate the allowed parameter range and to illustrate different aspects of the 
surviving points, we show the following plots (for positive sign of /i): 

1. In Fig. 1, the colour codes illustrate the percentage of Higgsino-like neutralino on 
the left plot and the neutralino contribution to the thermal DM relic density on 
the right. Note that in most of the points the neutralinos contribute only a small 
percentage of the total DM relic density and other DM candidates, such as axions, 
have to be present. The thermal DM relic density is close to the observed Planck 
value only in a small region corresponding to an LSP that is approximately 50% 
Bino and 50% Higgsino. However, as we will see in the next section, after imposing 
indirect detection constraints from Fermi, the only region which survives is the one 
where thermal DM is underabundant (by about 10% of the observed relic density). 
On the other hand, as we shall see in the next section, non-thermal Higgsino-like DM 
can lead to larger relic densities which can saturate the Planck value for reheating 
temperatures around 2-3 GeV. 

2. The colour codes in Fig. 2 illustrate the dependence on tan/3 on the left plot and 
different values of gaugino masses on the right plot with well-defined domains for 
different ranges of gaugino masses inside the V-shaped band. Note that tan/3 tends 
to have larger values as expected from the fit of rrih and the parameter ranges in this 
scan. Smaller values of A/M and vnjM are preferred for larger gaugino mass due to 
RG flow of masses to fit the experimental value of 

3. For Fig. 3 the colours illustrate on the left the different values of the typical scale 
of SUSY particles Msusy, defined here as the averaged stop mass Mg^gy = 

Notice that MgusY is around 4-5 TeV. In principle, we could explore values larger 
than 5 TeV however it would bring us beyond the level of applicability of the spectrum 
generator SPheno we have been using which assume similar values for all soft terms. 
An analysis for a split-like SUSY case with larger differences between sfermions and 
gaugino masses would be required in that case but this goes beyond the scope of this 
article. The colours on the right plot illustrate the dependence on the Higgs mass for 
which we have taken rrih — 125,125.5,126 GeV respectively. Note that for rrih = 126 
GeV there are allowed points only on the left of the V-shaped band because of the 
above mentioned cut-off on MgusY- Generally speaking we see that by allowing a 
larger range for the Higgs mass, we widen the V-shaped region. 

3.3 Direct and indirect detection constraints 

In the figures above we have set /i > 0 but their pattern does not change for /i < 0. The 
next step is to impose the following phenomenological constraints for the separate case of 
positive and negative /i since the DM direct detection cross section depends on sign(/i): 

• Fermi bounds on DM indirect detection [20]; 

• IceCube [23] and XENONIOO [21] bounds for spin dependent DM direct detection; 

• LUX [24], CDMS [22] bounds on spin independent DM direct detection. 
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Figure 3. Correlation between a = m/M and b = A/M for different values of the averaged stop 
mass Msusy (left) and the Higgs mass (right) in the region where the LSP is at least 10% Higgsino. 

3.3.1 Results for positive /i 

If we impose the above constraints, indirect detection bounds turn out to be very severe. 
Using the new Fermi bounds (pass 8 limit) coming from data collected until 2014 (the pass 
7 limit includes only data until 2012), we do not find any allowed point for Tr <2 GeV. 

In Fig. 4 we show the results for different reheating temperatures. The red points 
show the parameter space where we saturate the DM content measured from Planck [19]. 
We find more allowed points for Tr — 10 GeV compared to Tr = 2 GeV since the ratio 
of Tf/TR becomes smaller for large Tr, and so a smaller annihilation cross section is 
needed, resulting in a better chance to satisfy the bounds from Fermi. We only see the 
large Higgsino dominated regions for smaller Tr since in this case a larger annihilation 
cross section is needed to saturate the DM content. For larger Tr, an LSP with a smaller 
Higgsino component becomes allowed (higher Bino component) and this region appears for 
smaller values of gaugino mass. We finally mention that the Planck constraints on indirect 
detection through DM annihilation during the recombination epoch are less stringent than 
those of Fermi for our scenarios. 

Fig. 5 shows the spin independent and spin dependent WIMP-nucleon cross section 
after imposing LEP, LHC, Planck and Fermi constraints for Tr = 10 GeV. IceCube bounds 
rule out the orange and red regions of the right-hand side plot of Fig. 5. 

Fig. 6 shows the inclusion of LUX bounds on the spin independent direct detection 
constraints which rule out most of the points in Fig. 5 apart from a region corresponding 
to LSP masses around 300 GeV which is at the border of detectability. 

Moreover, for Tr = 2 GeV there are red points which saturate the observed DM 
content. In this case the neutralinos are becoming more pure Higgsinos in order to enhance 
the annihilation cross section and the Fermi constraint becomes harder to avoid but there 
are still regions allowed by both direct and indirect detection searches. There are more 
green points for Tr = 10 GeV since the annihilation cross section becomes smaller due to 


- 9 - 










5 10 



m^= 125.5- 126 GeV 
Tf,= 2GeV 




I I I I I I I I 


5 

A/M 


S 10 




mrt= 125.5-126 GeV 
Tn= 2 GeV 



• •.*••. 

"•Vi£,k&»r : 



= 125.5 - 126 GeV 
Tp = 10 GeV 




5 

A/M 



Figure 4. Case with /i > 0: m/M vs A/M after imposing LEP, LHC, Planck and Fermi constraints 
(left) and corresponding LSP composition (right) for rrih = 125.5 - 126 GeV and Tr = 2,10 GeV. 


Bino mixing which means a larger allowed region after using Fermi data but the constraint 
from the direct detection becomes more stringent (due to Bino-Higgsino mixing in the LSP) 
and so there are no red points which saturate the observed DM content.^ For the points 
shown in Fig. 6, the GUT values c — B/M and d = ji/M are around 0.6 and 1 respectively. 

3.3.2 Results for negative /i 

The results for the negative fi case are shown in Figs. 7 and 8. Note that under the same 
conditions the allowed parameter space for /i < 0 is larger than the one for /i > 0. This 
can be understood as follows [37]: o-^o_p is dominated by the t-channel /i, H exchange 
diagrams which mostly arise from down type (s-quark) interaction: 


oc rrid 


cos a Fh sin a 


cos j3 m? 


H 


cos P ml 


(3.8) 


®The direct detection exclusion however depends on various uncertainties, e.g. strange quark content of 
proton, form factor etc. [38]. 
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Figure 5. Case with /i > 0: correlation between a — m/M and h = A/M after imposing LEP, 
LHC, Planck and Fermi (pass 8 limit) constraints with the corresponding spin independent (left) 
and spin dependent (right) WIMP-nucleon cross section for rrih = 125.5 GeV and Tr = 10 GeV. 
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Figure 6. Gase with > 0: correlation between a = m/M and h = A/M after imposing LEP, 
LHC, Planck, Fermi and LUX bounds for = 125.5 - 126 GeV and Tr = 2,10 GeV. 


where a is the Higgs mixing angle, = (A^i 2 — A^n tan 0vp)(^i4 cos cr + A^issincr) and 
Fh = (A^i 2 -A^ii tan 0 vp)(A^i 4 sincr-A^is coscr) using x? = iViiH + A^i2lU+A^i3^i + A^i4^2- 
For /i < 0, the ratio A^i 4 /A^i 3 is positive and this amplitude can become small due to 
cancellations if: 

A^14 tan a + m^/cot a 


(3.9) 


1 + ml/m% 

is satisfied (for tana < 0). In Fig. 8 there are more allowed points compared to Fig. 6 even 
if there are still no points which saturate the observed DM content for Tr = 10 GeV due 
to stringent direct detection constraints. For the points shown in Fig. 8, the GUT scale 
values of B and fi are still both of order M. 
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Figure 7. Case with /i < 0: a versus b after imposing LEP, LHC, Planck and Fermi (pass 8 limit) 
bounds (left) and LSP composition (right) for rrih = 125.5 - 126 GeV and Tr = 2,10 GeV. 


*: j* iJ. I * 


125.5- 126 GeV 
Tf,= 2GeV 




llllllllllllllllllllLllIlllLillll 


-4-3-2-101234 


5 

A/M 



Figure 8. Case with /i < 0: correlation between a — m/M and h = A/M after imposing LEP, 
LHC, Planck, Fermi (pass 8 limit) and LUX for rrih = 125.5 - 126 GeV and Tr = 2,10 GeV. 
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4 Discussion of results 


4.1 Analysis of the allowed parameter space 

We are now in a position to put all our results together and explain the effect of each ex¬ 
perimental bound on our parameter space. In the end we shall analyse the spectrum of su¬ 
perpartners that appears for the points surviving all the phenomenological constraints. All 
the observables analysed in this section have been computed numerically using mi cr OMEGAS 
v3 [48] linked to SPheno v3.3.3 [47]. 

Fig. 9 shows the relation between the spin independent WIMP-nucleon cross section 
and the LSP mass. Depending on Tr there is a different upper bound for neutralino masses 
which is given by the Planck constraint on DM. For larger values of Tr, the non-thermal 
relic density decreases, and so heavier neutralinos can pass the Planck constraint on the 
DM relic density. On the other hand, for larger Tr the parameter space for standard 
thermal DM (orange band) becomes also larger. Note that LUX 2013 results exclude at 
90% most of the parameter space and the next round of results (LUX 300 days) will be 
able to probe the remaining regions (the light blue points below the LUX line). 

In the scenario we considered, the gaugino masses are unified at the GUT scale and 
therefore the evolution of electroweakinos is totally dominated by the RG ffow. This 
implies that the LSP can only be Higgsino- or Bino-like (or a mixed combination of them). 
The largest contributions to the thermal averaged annihilation rates are given by (see for 
example [15] and references therein): 

(<Teffv) = + 3tan^6»i4/ + Ow) , (4.1) 

for Higgsino-like neutralinos (in the limit Mw "C /i) annihilating into vector bosons through 
chargino or neutralino interchange, and: 

g2^&v?9w -Q/)^r(l + r^) 

^ 27vm?~{l + r)^ 

for Bino-like LSP annihilation into fermion-antifermion (Ta^ and Qf are the third compo¬ 
nent of isospin and the fermion charge and r = Ml/m?-). This process is driven at tree 

level by the t-channel exchange of a slepton /. In the case where the LSP is a mixed com¬ 
position of Higgsino and Bino, the expression of the annihilation rate is an interpolation 
between (4.1) and (4.2). Fig. 9 shows also the effect of Fermi bounds. As suggested by 
(4.1) and (4.2), the most constrained regions are those with smaller LSP masses. The grey 
band corresponds to points excluded by LEP bounds on chargino direct production. 

Fig. 10 shows the amount of non-thermal DM relic density provided by the LSP in 
terms of its mass, together with the bounds from indirect detection and LUX. The Planck 
value of the DM content can be saturated in the region which is not ruled out by direct 
detection bounds only for Tr — 2 GeV. Given that for larger Tr the amount of LSP 
DM gets smaller, the cases with Tr > 2 GeV require multi-component DM. Combining 
Fig. 9 and Fig. 10, we find that the LUX allowed regions, indirect detection limits and 
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the abundance of LSP DM are correlated for different Tr. The allowed regions, where the 
observed DM content is saturated, depend on Tr but they are generically around ^ 300 
GeV. 
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Figure 9. Spin independent cross section versus LSP mass for /i > 0 (left) and /i < 0 (right). The 
light blue points are not ruled out by indirect detection experiments. We show up to 800 GeV. 
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Figure 10. Non-thermal DM abundance predictions versus LSP mass for /i > 0 (left) and /i < 
0 (right) and Tr = 2 , 5,10 GeV. The light blue points are not ruled out by indirect detection 
experiments while only the yellow points are allowed by LUX 2013 results. 
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In Fig. 11, we show the Planck constraints on indirect detection through DM anni¬ 
hilation during the recombination epoch. We use WW final states corresponding to an 
efficiency factor /eff = 0.2 [39]. Even going all the way to the cosmic variance bound, these 
constraints turn out to be less stringent than those coming from Fermi and LUX. On the 
other hand. Fig. 11 shows the correlation between the LSP composition and the bounds 
coming from both direct and indirect detection. Concerning Fermi and Planck limits on DM 
annihilation, these bounds allow almost all possible combinations of Bino/Higgsino neu- 
tralinos. The restrictions coming from the ff and WW channels depend on the neutralino 
composition: for a Higgsino-like LSP, the most stringent constraint is due to annihilation 
into vector bosons, while for a Bino-like LSP the main constraint comes from annihilation 
into a fermion-antifermion pair. 

The LUX constraints in Fig. 11 reduce the parameter space to the region where the 
LSP is mostly Higgsino-like. This could be a bit puzzling since the WIMP-nucleon cross 
section is dominated by the Higgs exchange channel: 


- ff 

^x-P 0 ^- -4 - 


(4.3) 


where a^, and are respectively the Higgsino, Bino and Wino LSP components. 
According to this expression, the cross section is enhanced when the Higgsino component 
increases. However in Fig. 11 direct detection bounds allow only points which are mainly 
Higgsino-like. The reason of this effect is in the effective coupling yy/i which for a Bino-like 
LSP looks like: 


^XXh — 


mz sin 9w tan 9w 
Ml - /i2 


(Ml -h /i sin 2/3), 


(4.4) 


where for moderate to large tan/3 the second term is negligible and ji > Mi. Hence this 
coupling is dominated by Mi. On the other hand, the coupling for a Higgsino-like LSP is: 


Cyy/i ^ - {1± sm2P) (tan^ 9w 


mz cos 9 ^ mz cos 9 

Ml - l/x A M2 - ImI 


(4.5) 


where ± is for the and components and ji < Mi. Contrary to the Bino-like 
case, this coupling is inversely proportional to Mi. Thus the WIMP-nucleon cross section 
grows in the regions where the LSP is Higgsino-like and Mi is small or where the LSP 
is Bino-like and the gaugino mass is large. If we compare Fig. 2 (right) which shows 
the distribution of gaugino masses along the V-shaped band, with Fig. 4 (right) and 7 
(right), we realise that the region with a smaller cross section is the one at the bottom 
of the V-shaped band where gaugino masses are big and the LSP is very Higgsino-like. 
The region where the LSP is more Bino-like has smaller gaugino masses and the cross 
section is larger. Fermi constraints however become more stringent in the case with more 
Higgsino content due to larger annihilation cross section. The competition between LUX 
and Fermi constraints produces the allowed parameter space where the Planck value of the 
DM content is saturated for Tr = 2 GeV. 
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Figure 11. Comparison between detection constraints from Planck, Fermi and LUX for /i > 0 
(left) and /i < 0 (right). We have set /eff = 0.2. We show up to 900 GeV. 

Finally, in Fig. 11 there is a change of behaviour of the thermal averaged cross section 
for masses around 130 GeV. The reason is the following: this region is (as can be shown in 
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the plot) Higgsino-like, but the masses are closer to M\y and {av) is no longer described 
by (4.1) but by something like (with x = jj^lmw)- 


In Fig. 12, we show the spectra of SUSY particles for the allowed regions of Fig. 6 (blue 
points below the LUX line). We find that sleptons, staus, Higgses, all other scalar masses 
and gluinos are rather heavy since they are between about 2 and 7 TeV. The lightest and 
second to lightest neutralino and the lightest chargino are around 280 - 340 GeV while all 
other neutralinos and charginos are heavy. The allowed region for Tr = 2 GeV is shown 
on the left side of the vertical line with the label Tr — 2 GeV where the points situated 
exactly on the line satisfy all the constraints including the current DM content as measured 
by Planck. Similarly, the allowed region for Tr > 5 GeV is shown on the left side of the 
vertical line with the label Tr > 5 GeV even if there are no points in this region which 
saturate the current DM content. Notice that the spectrum is essentially independent of 
the reheating temperature Tr and the hierarchy between the different sparticles is robust. 



Figure 12. The mass spectra of superpartners for allowed points shown in Fig. 6 for different 
values of Tr. 


4.2 Astrophysical uncertainties 

The direct detection cross section can involve various uncertainties, e.g. strange quark 
content of proton, form factor, local DM density and LSP contribution to the total amount 
of observed DM abundance. The local density can be 0.1-0.7 GeV/cm^ [40]. There could 
also be astrophysical uncertainties in the indirect detection results beyond what has been 
considered so far. Recently, it is mentioned in [41] that if the thermal neutralinos do not 
produce the entire amount of cold dark matter, the direct and indirect detection cross 
sections should be reduced by R and respectively with R = D/i^/0.12. Possible bounds 
arising from Fermi are now almost negligible since they are suppressed by R?. Once the 
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suppression factor R is taken into account, Fermi, Planck and other indirect detection 
experiments have lower impacts. Concerning the effect on LUX and other direct detection 
bounds, the cross section is now reduced by R which is equivalent to multiplying the 
effective couplings (4.4) and (4.5) by \fR. This clearly introduces a new parameter in the 
discussion performed in the previous section. 

If we assume such a reduction in the cases where < 0.12, more parameter space 

could be allowed for multi-component DM regions as shown in Fig. 13. The pink region 
is disallowed by Fermi data. The Tn dependence of the region constrained by Fermi in 
Fig. 13 is due to the fact that the factor R is now a function of Tr 

^ Tf 

R becomes larger for smaller values of Tr {—2 GeV) and the Fermi constraint becomes 
important. 

The region below the dashed line satisfied present LUX limits. In particular this implies 
that the region with lighter neutralinos is now unconstrained by LUX. This region typically 
corresponds to more Bino component in the LSP as shown in Fig. 4 (right) and 7 (right). 
We have therefore a different situation compared to before, because now neutralinos with 
a larger Bino component are allowed. 

Let us stress, however, that the prediction for the region where saturates the 

DM content remains unchanged, i.e. only the case Tr — 2 GeV contains points which are 
still allowed by all data and saturate the DM content with an LSP mass around 300 GeV. 

This new factor R helps us to extract more parameter space for the multi-component 
DM scenarios. However, the DM simulations need to establish the validity of the assump¬ 
tion that proportions of various DM components in the early universe is maintained even 
after the large scale structures are formed. 
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Figure 13. R x _~o vs m^o for Tr = 2, 5 GeV. R is defined in the text. More parameter space 
is allowed. 
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5 Conclusions 


Non-thermal DM scenarios emerge in UV theories like string theory due to the presence 
of gravitationally coupled scalars which decay at late times when they are dominating 
the energy density of the universe. In such models the reheating temperature due to 
moduli decays is typically below the freeze-out temperature, Tr < Tf when assuming 
an MSSM particle as main DM constituent. In this paper, we have studied the non- 
thermal version of the CMSSM/mSUGRA and contrasted it with both particle physics 
and astrophysical experimental constraints. The experimental information available at 
this moment, including the well known value of the Higgs mass, is enough to rule out the 
vast majority of the non-thermal CMSSM parameter space. There is still a small region 
which is consistent with all observations and is at the edge of detection by both astrophysics 
and particle physics experiments, resulting in a very interesting situation for beyond the 
SM physics. 

In our determination of the allowed parameter space we have used constraints from 
collider experiments (LEP, LHC), cosmic microwave background observations (Planck) and 
direct and indirect DM searches (Fermi, XENONIOO, IceCube, LUX, CDMS). We have 
found that the most restrictive constraints come from Fermi and LUX which single out 
a small region of the non-thermal CMSSM parameter space corresponding to a Higgsino- 
like LSP with a mass around 300 GeV that can saturate the observed DM abundance 
for Tr ^ 2 GeV while larger reheating temperatures require additional contributions to 
the present DM abundance. These results are summarised in Fig. 10 which shows the 
comparison between the cases of positive and negative /i. 

This non-thermal scenario leads to a clear pattern of SUSY particles. In particular, 
the fact that the LSP is Higgsino-like makes the lightest chargino, the lightest neutralino 
and next to lightest neutralino to be almost degenerate in mass. This kind of scenario can 
be probed at the LHC using monojet plus soft leptons plus missing energy [42], monojet 
signal [43] and two Vector Boson Fusion jets and large missing transverse energy [44]. On 
the other hand, all the other superpartners are much heavier and beyond the LHC reach 
but accessible to potential future experiments such as a 100 TeV machine. 

It is worth mentioning that non-thermal CMSSM scenarios with TeV-scale soft terms 
and reheating temperatures around 1-10 GeV can emerge in string models where the 
visible sector is sequestered from the sources of SUSY breaking [11, 12]. In a subset of 
the parameter range these string scenarios lead to Mgoft ^ Mpe^ <C mmod <C 

^GUT Mp^!^ and Tr ^ Mp^!^ where e m^i 2 lMp <C 1. For e ^ 10“^, one obtains 
TeV-scale soft terms, Mqut ^ 10^^ GeV, Tr ^ 1 — 10 GeV and mmod 10^ GeV for 
ms /2 - 10^^ GeV. 

We point out that our analysis is based on the CMSSM/mSUGRA for which all su¬ 
perpartners are expected to be at similar masses close to the TeV scale. In this sense we 
restricted ourselves to scalar masses lighter than 5 TeV which is the range of validity of the 
codes we have used to perform our analysis. There are however several ways to generalise 
this simplest scenario: 

• Consider non-universal extensions of the CMSSM. Small departures from universality. 
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even though strongly constrained by flavour changing neutral currents, allow more 
flexibility in the parameter space and will slightly enhance the allowed region. 

• Consider sfermions heavier than 5 TeV as in the split SUSY case. This is not only an 
interesting phenomenological possibility but is also the other class of scenarios that 
were derived in the string compactiflcations studied in [11, 12]. 

• Consider MSSM scenarios with R-parity violation and late decaying moduli flelds. 
This class of models has not been studied in detail and given the fact that even with 
R-parity conservation there seems to be a need for other sources of DM such as axions 
or axinos, this should be a possibility worth studying. 

• Consider explicit D-brane models which tend to generate models beyond the standard 
MSSM (see for instance [45, 46] for such models). 

It is encouraging that new planned experiments such as upcoming LUX result and 
XENONIT will be enough to rule out the rest of the allowed parameter space, independent 
of the upcoming LHC run. Clearly also new LHC runs and future planned colliders will 
be crucial for this class of models. Combining astrophysical and collider measurements is 
probably the most efficient way to constrain beyond the SM physics and this article is a 
clear illustration of this strategy. 
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